Pairing on striped t-t'-J lattices 
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Results are given from a density matrix renormalization group study of pairing on a striped t- 
t'-J lattice in the presence of boundary magnetic and pair fields. We find that pairing on a stripe 
depends sensitively on both J/t and t' /t. In the strong-pairing model-parameter regime the stripes 
are easily coupled by the pair field, and have a uniform phase. There is a small but measurable 
energy cost to create anti-phase superconducting domain walls. 
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Experimental studies are providing new insight into 
the interplay of the charge, spin and d-wave pairing corre- 
lations in the underdoped cuprates. Scanning tunneling 
microscopy measurements Q on Cai.gsNa 0.12CUO2CI2 
and Bi2Sr2Dyo.2Cao.8Cu208+8, suggest that at low tem- 
peratures, as the doping increases, superconducting cor- 
relations develop on a glassy array of 4ao wide do- 
mains oriented along the Cu-0 x or y-bond direc- 
tions. Neutronfi!] and x-ray scattering[3| experiments on 
Lai.873Bao.i25Cu04 find charged stripes with a 4ao pe- 
riod separated by 7r-phase shifted antiferromagnetic re- 
gions. Moreover, when the temperature decreases be- 
low the spin ordering temperature, the planar Pnb re- 
sistivity shows evidence of a Kosterlitz-Thouless[J| like 
behavior consistent with the development of 2D pair- 
ing correlations.^] Remarkably, Pab follows the Halperin- 
NelsonQ 2D prediction over an extended temperature 
range, implying a decoupling of the pair phase between 
the Cu02 planes. In underdoped La2-xSr2;Cu04, super- 
conductivity and static spin density waves coexist 0], and 
recent far-infrared measurements [8] find that the Joseph- 
son plasma resonance is quenched by a modest magnetic 
field applied parallel to the c-axis. An applied c-axis 
magnetic field is known to stabilize a magnetically or- 
dered state for a range of dopings near x = 1/8. This 
is believed to be a striped state and Schafgans et aZ.Q 
have argued that just as in Lai 873Bao.i25Cu04, the es- 
tablishment of antiferromagnetic stripes leads to a sup- 
pression of the interlayer Josephson coupling. To explain 
this suppression, it has been suggested that anti-phase 
domain walls in the d-wave order paramater, locked to 
the SDW stripes, are stablized in the striped magnetic 
state. [9I, Kd\ In this case, the 90° rotation of the stripe 
order between adjacent planes would lead to a cancella- 
tion of the interlayer Josephson coupling. 



Early mean field calculations [llj found striped states 
in t-J and Hubbard models. However, these stripes 
had a hole filling which was twice that which was ob- 
served in the cuprates. While arguments were made 
that this problem could be overcome by including a next 



near neighbor hopping t' , an alternative view suggested 
that it reflected the importance of underlying dx2_y2 
pair field correlations. 13] Density matrix renormaliza- 
tion group (DMRG) calculations found half-filled hole 
stripes, 7r-phase shifted antiferroniagnetism and short 
ranged d^2_y2 pairing correlations. 14] This interplay of 
oscillating hole density, spin density and (1^2 _y2 super- 
conductivity was also found in Gutzwiller- pro jected vari- 
ational Monte Carlo (VMC) calculations. [lOj In certain 
parameter ranges, these VMC calculations found a stable 
striped state in which the (1^2 _y2 pair field had tt phase 
shifts between the stripes, i.e. anti-phase domain walls. 
A recent renormalized mean field theory (RMFT) treat- 
ment found a similar small energy difference, with the 
uniform phase d-wave state lying slightly lower in energy. 
jl^J A TT- phase shifted d-wave pair field would provide a 
natural explanation for the observed suppression of the 
interlayer Josephson coupling. (9l.[lo| A similar, low lying, 
modulated superconducting state was also found in VMC 
resonating valence bond (RVB) calculations. 3 121 Here, 
however, no incommensurate AF order was assumed. fl3| 
We have carried out a series of DMRG calculations on 
under-doped t-t'-J lattices with a Hamiltonian 

H = -i +h-c.) - t' T.i^j)MsCjs + h.c.) 
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Doubly occupied sites are excluded from the Hilbert 
space, Si and c| ^ are electron spin and creation oper- 
ators, respectively, and Ui = c|^|Cif +cl^Cii is the electron 
number on site i. There is a near neighbor (ij) hopping 
t, a next near neighbor (ij)' hopping t' , and an exchange 
coupling J. We set t = 1. Using boundary magnetic and 
pair fields, we have explored how a pair field is estab- 
lished over a magnetically striped array at low doping. 

Fig. [T] shows two 12 X 8 ladders (tubes) with cylindrical 
boundary conditions (CBCs). In Fig. [Ija), a weak stag- 
gered magnetic field is applied to the open ends and the 
number of holes is fixed at 12, corresponding to a doping 
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FIG. 1: (a) Hole (1 — rii) and spin densities (S'f ) for a 12 x 8 
lattice with 12 holes and J = 0.5 and t' — 0, with cylindrical 
boundary conditions: periodic in the y-direction, open in the 
a;-direction. A staggered magnetic field of magnitude h = 
±0.05 has been applied to the ends (red X's). (b) A similar 
lattice with = 0.55, Jy = 0.45 and t' — 0.0 such that the 
stripes run along the a;-axis. Here we use h = ±0.2 on the 
sites with a red X, and a pair field Ao = 1.0 to the edge links 
without X's. A chemical potential jj, = 1.23 was used to give 
a doping of x = 0.127. (c) The pair field strength (-Dij) on 
each link for the system shown in (b). 



d = 0.125. The end boundary conditions break the trans- 
lational and spin symmetries in the x-direction, giving 
rise to finite, varying values of (n^) and (Sf), but the ba- 
sic stripe pattern, with antiferromagnetic order between 
stripes and with tt phase shifts in the magnetic order 
across stripes, is intrinsic and the open ends and bound- 
ary h field only act to pin the stripes. Up to m = 4000 
states per block were kept. 

In this cluster there are four holes per stripe, and 
DMRG calculations on longer stripes have shown that 
this linear stripe filling of 0.5 holes per unit length is the 
preferred filling. 1J| Fluctuations in which a pair of holes 
are exchanged between the stripes are energetically unfa- 
vorable for this short a stripe, and the local dx2-y2 pair 
field correlations are short ranged, and pair field coupling 
between stripes is negligible. 

Previous DMRG calculations have not found ground 
states with both extended pairing correlations and 
stripes. There are two reasons: first, it has been diffi- 



cult to construct limited-size clusters allowing significant 
particle number fluctuations on a stripe, and second, as 
we show below, the model parameters which strongly fa- 
vor pairing (e.g. J/t ~ 0.5, t' /t ~ 0.2) are different from 
the values usually taken to represent the cuprates (e.g. 
J/t ~ 0.3, t'/t = -0.2). In Fig. [Hb), we show results 
for a cluster which does have both stripes and pairing. 
In order to allow hole fluctuations, a slightly anisotropic 
exchange interaction {J^ — 0.55, Jy = 0.45) was cho- 
sen to favor orienting the stripes along the x-direction, 
overcoming an opposite tendency due to the cylindrical 
geometry. Then, in addition to the magnetic fields at the 
open left and right ends, a pair field coupling has been 
applied to the ends of the stripes. Defining the link pair 
creation operator 
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on specified links, we add to the Hamiltonian boundary 
region terms of the form Ao-Dij, where 
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and we measure Dij on each link in the resulting (approx- 
imate) ground state. For this system we took Aq — 1.0 
for the four thickest lines of (c), and also Aq = 0.5 for 
the four vertical links adjacent to them. With the pair 
field boundary conditions, total particle number is only 
conserved modulo two, and the average number of holes 
is controlled by a chemical potential fi. Up to to = 6000 
states were kept per block during 23 sweeps in this cal- 
culation. As shown in Fig.[ljc), a proximity d-wave pair 
field is established throughout the lattice. 

In order to understand the system in more detail, it is 
useful to separate questions dealing with (a) pairing on 
a stripe from (b) pairing between stripes. First, can a 
single stripe support strong pairing, and if so, for what 
model parameters? Second, do stripes with pairing cou- 
ple their pair fields, and if so, is the coupling in phase 
or antiphase? To answer these questions we will study 
clusters somewhat smaller than shown in Fig. 1, to ease 
the computational burden and increase the accuracy. 

To look in more detail at the pairing correlations as- 
sociated with a stripe, we have studied a single stripe on 
the 16 X 5 lattice with CBCs, shown in Fig. In this 
case boundary conditions were used to force the presence 
of a stripe similar to those shown in Fig.[T](b)-(c). Here a 
strong pair field was applied to only one end of the stripe. 
The proximity induced pairing response (Dij) is shown in 
Fig. IHb) for a case with strong pairing. The strength of 
the induced pair field depends upon J, t' and the doping 
X. As a measure of this response, its magnitude on the 
12th y = 2 — 3 rung is plotted in Fig.[2I^c) versus the hole 
doping for various values of t' and J. As previously found 
in both DMRG and VMC calculations, a positive value 
of t' favors pairing while a negative value suppresses it. 
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FIG. 2: A cluster with a single forced stripe on a 16 x 5 
ladder with CBCs, with the x direction oriented vertically. 
For all dopings, on the y = 5 leg, a staggered field { — l)^h 
with h — 0.05 was applied, along with a chemical potential of 
2.0. On the other 4 legs, a chemical potential /i was applied 
to vary the doping; in the case shown in (a) and (b), fi — 1.41, 
J — 0.5, and t' — 0.2, yielding a doping of x = 0.106, which 
corresponds to a linear doping of 0.53. In each case a strong 
pair field was applied to four rungs {x — 1 — 4, connecting 
y = 2 and y = 3), visible as the four thickest links in (b). 
The magnitude of the applied field in the four rungs was 1.0, 
1.0, 0.5, 0.25. (a) The hole density (1 — rii) and spin density 
(Si), (b) The measured pair field (Dij). (c) The measured 
pair field at the a; = 12, y = 2 — 3 rung versus the number of 
holes per unit length. The maximum value for {Dij) in (c) is 
of order half of what one would find for a BCS ground state 
with Ao/t = 0.1. 



Here one also sees that when the pairing is strongest, the 
response peaks for a linear filling p ~ 0.5 holes per unit 
length, but shifts to higher doping for smaller t' . We see a 
strong dependence on J/i, with pairing when J/t = 0.3 
quite weak and when J/t = 0.5 quite strong. Also, as 
shown in Fig. [31 the compressibility, which is related to 
the slope of the curves, for p = 0.5 increases as t' in- 
creases, consistent with the observed enhancement of the 
pairing response for positive values of t' . The increased 
pairing for larger J/t may be due to a reduced repulsion 
between pairs, leading to an enhanced compressibility. 

Next we turn to the question of anti-phase domain 
walls in the pair field. For values of t' and doping where 
there is a significant pair field response, we find that the 
pair field remains in phase across the stripes, as shown 
in Fig. [TJb). This implies that the energy to create an 
anti-phase domain wall is positive. To probe this, we 
have studied 12 x 6 ladders with open boundary condi- 
tions in both the x and y directions, and with magnetic 
and chemical potential fields applied to force two stripes. 
Fig. [4] shows such a system where each stripe has a linear 
doping of 0.6. Then by applying pair fields on every link 
on the outermost legs {y — 1 and y — Q), a pair field is 
established. When the applied pair fields are in phase, 
the induced pair field is shown in the middle figure of 




FIG. 3: The linear density versus p for different values of the 
next nearest neighbor hopping t' for J = 0.5 for the systems 
of Fig. 2. 




FIG. 4; (a) Hole and spin densities for a 12 x 6 ladder (plotted 
with X oriented vertically) with J = 0.5 and t' = 0, and with 
open boundary conditions in both directions with fields ap- 
plied to force two stripes. The center two legs have an applied 
staggered field ( — l)^/i with h = 0.05, along with a chemical 
potential of 1.3. The outer four legs have a chemical potential 
of 0.9, leading to a doping of x=0.1981, corresponding to a 
linear doping per stripe of 0.594. The outermost two legs have 
applied pair fields of 0.5 on each horizontal rung. In (a), the 
pair fields are applied in phase, (b) Measured pair fields for 
the in-phase case, (c) Measured pair fields when the applied 
fields have opposite sign. 



Fig. [4] The right hand figure shows what happens when 
the applied pair fields are out of phase. The energy per 
unit length versus the DMRG sweep number for the in- 
phase and anti-phase configurations are plotted in Fig. O 
for different values of t' . Although the energy difference 
is small, it varies little with the sweep as the number of 
states in increased, up to to = 3000 for sweep 17. Here 
one sees that it costs energy to create an anti-phase do- 
main wall and the energy per unit length increases as 
the overall strength of the induced pair field increases. 
For t' = 0.2, the energy per unit length of the anti-phase 
domain wall is of order O.Olt. 

There arc both similarities and differences between 
our DMRG results and those from VMC.0 Both ap- 
proaches find evidence for low lying striped states with 
d^2_y2 pair fields, as seen experimentally. However, we 
find that negative values of t' suppress the c?-wave pair 
field. Thus the parameter regime where we have studied 
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FIG. 5: Energy per site for the 12 x 6 ladder of Fig. 4 versus 
sweep number for the in phase and antiphase applied fields. 
The runs all had J = 0.5 and had the values of t' indicated. 



the interplay of the stripes and the pair field differs from 
the t' < region of the VMC study. Our DMRG cal- 
culations also find that the energy to form an anti-phase 
d-wave domain wall is positive. 

In comparing to experiments for the cuprates, the puz- 
zle raised by our calculations, as well as the RMFt[i51| 
results, is the positive energy required to create an anti- 
phase domain wall in the d-wave order. Because this is a 
small energy, one can imagine that effects missing from 
the t-t'-J model could lead to antiphase domain walls. 
However, we believe that there could be an alternate ex- 
planation of the suppression of the interlayer Josephson 
coupling observed in the underdoped regime. It could be 
that the decoupling arises from the lack of overlap of the 
Fermi surfaces of the adjacent layers. There are a num- 
ber of experiments 3, 20] which imply that a Fermi sur- 
face reconstruction occurs for hole doping near x — 1/8. 
The resulting Fermi surface is characterized by electron 
pockets and open orbits whose location in the Brillouin 
zone depends upon the stripe orientation. [2l| In this case, 
since the stripes alternate their orientation by 90° from 
one plane to the next, the lack of overlap between the 
Fermi surfaces can lead to a suppression of the interlayer 
pair transfer processes (2^. 
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